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Abstract 


A quantum version of the Lighthill equation that originated in the 
field of theoretical aeroacoustics is derived for the probability density of 
a superfluid starting from the time-dependent Gross-Pitaevskii equation. 
It involves a second-order time derivative and should be supplemented by 
two-time boundary conditions. Its physical implications are discussed in 
relation to the quantum equilibrium hypothesis and the general applica- 
bility of Born’s rule. 


The hydrodynamic formulation of quantum theory has a long history, starting 
with Madelung [1| almost immediately after the discovery of quantum mechan- 
ics and coming to maturity with the work of Takabayasi and others [2] [3] [4] [5] 
6} [71 [8] [9]. According to this picture there is a close analogy between motion of 
material and probability fluids which is made apparent by two continuity equa- 
tions of the Navier-Stokes type for the probability density and the momentum 
of this probability field. This indicates that quantum phenomena may be mod- 
eled in a hydrodynamic framework. It should be remembered on the other hand 
that the material fluids exist in the three dimensional space but a many particle 
wavefunction in an abstract configuration space. Furthermore, the stress tensor 
in the case of a probability fluid has distinctive quantum features. This alter- 
native formulation of quantum mechanics, which shares many similarities with 
the pilot wave theory [8], did not prove especially popular among physicists. 
Yet the hydrodynamical picture is particular successful when it comes to the 
description of ultra cold atoms and weak interacting quantum gases [IOLI] even 
though the physical content is rather different in that occasion, since the macro- 
scopic single-particle complex wavefunction that describes a condensate is not 
a quantum state and the equation of motion is non-linear. The time-dependent 
Gross-Pitaevskii equation that governs the time evolution of a quantum super- 
fluid is written as 
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where g is a coupling constant, V(r,t) the external potential and W a single- 
particle wavefunction. 

A brief digression is befitting here in order to provide some motivation for 
what will follow. An initial observation is that the non-linear wave equation (1) 
contains a first order derivative in time so it resembles a diffusion equation in 
that respect which means that W(t) Æ W(—t). According to the Wigner pre- 
scription U(t) — W*(—t) it follows that the particle density is time symmetric 
n(t) = n(—t). It seems reasonable then to investigate if we could derive a man- 
ifestly time-invariant evolution equation for n(t), possibly second order in time. 
In other words we exploit the Gross-Pitaevskii equation for the macroscopic 
wavefunction to deduce an equation of motion for the probability density. To 
achieve this we proceed by introducing a Madelung polar transformation 


Y = yne's/h, (2) 


where S is the multivalued phase of the wavefunction and the velocity field 
coupled to a vector potential is expressed as 


= (3) 
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After separating real and imaginary parts two formulas are obtained. By equat- 

ing the imaginary parts we arrive at a probability density continuity equation 

which is written in index notation, where summation over repeated indices is 

implied, as 
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The Hamilton-Jacobi equation is obtained from (1) through (2) by equating its 
real parts. 
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In the above expression Q = -57 Tn ÍS the Bohm quantum potential. Taking 


the divergence gives the second hydrodynamic equation 


o 
dv) +V- (nvv)=-V(V +Q+gn), (6) 
where the quantum force is the divergence of the sum of both classical and 


quantum potentials. It is possible to write the above in a tensor form as in [10] 
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The probability quantum stress tensor is written as 
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Notice that the stress tensor does not depend on the velocity field but only 
on the density of the probability fluid. It is a second rank tensor with nine- 
components for a single particle in the three dimensional space, whereas for N 
particles the number of components is (N?. When g = 0 we recover the fa- 
miliar hydrodynamic equations that follow from the single particle Schrödinger 
equation [4]. Since we have performed differentiation the hydrodynami- 
cal equations we obtained are not equivalent to the Gross-Pitaevskii unless we 
impose certain constraints on the line integral of the velocity field (see the rel- 
evant discussion in [6] for the Schrédinger equation). As we do not aim 
to discard the wave equation or to recover the wavefunction equation from the 
hydrodynamic ones, this does not affect our scope. 

It is possible to rearrange and merge (4) and (7) into a single equation in 
a fashion similar to Lighthill’s who long ago derived an evolution equation for 
the generation of aerodynamic sound [15]. Differentiating over time 
(4) and taking the divergence of (7) and then subtracting the two resulting 
expressions we manage to eliminate the terms associated with the momentum 
of the probability fluid. To make more explicit the connection to the Lighthill 
equation we further subtract from both sides a c2V?n, where co an arbitrary 
reference speed, and arrive at the following time-invariant formula 
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where the expression inside brackets in the right hand side is a Lighthill source 
tensor which includes a quantum component with highly non-classical proper- 
ties. Hence we have obtained the equation we sought for n(t). This is our 
central result, an exact, non-linear tensor equation for the probability density 
that is driven by the velocity field governed by (1). We notice that the particle 
density appears also in the right hand side in a non-linear manner through the 
quantum stress tensor and other terms. Formally the solution can be written 
in integral form as 


n(x, t) By, (10) 
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where r = |x—y| the distance between reception and source and T;; the Lighthill 
tensor mentioned earlier. Notice, that we have included only the retarded solu- 
tion but advanced ones may be included too when dealing with appropriately 
pre- and post-selected quantum ensembles (in that case we require two-time 
boundary conditions for the particle density and its first time derivative). 
Due to the non-linear character of the source tensor in (9), even when g = 0 


which means that the driving equation (1) is linear, makes this equation difficult 
to solve for practical purposes yet it is proposed here for the conceptual insight 
it may offer to the study of foundational quantum mechanical theory. Some 
remarks should be included for the case with g = 0. Firstly, a solution Ų mul- 
tiplied by a complex number w yields a new solution wW. This illustrates the 
homogeneity property which is reflected into (9) when we multiply a solution 
n by a factor |w|?. Secondly, since the evolution equation for n is non-linear, 
a superposition of number densities does not constitute a new solution. Yet, 
a superposition of wavefunction solutions satisfies both the linear Schrödinger 
equation and (9) provided that g = 0. 

An equation second order in time derived from linearized hydrodynamic 
equations is given in [7]. We instead linearize (9) by letting n + no + ôn 
(no = |W/? is an equilibrium solution as will be explained in the final paragraph 
and ôn a small perturbation). Keeping first order perturbation terms yields 
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with a stress tensor perturbation 


2 
oli; = 5 [Vi V; = VilnnoVj _ VjlnnoVi + Vj lnnoVilnno] ôn. (12) 
A word of caution is appropriate at this point. The Gross-Pitaevskii equa- 
tion (or its linear counterpart when g = 0) is not to be replaced by (9). Neither 
do we seek to derive the Schrödinger equation (linear or non-linear) from it. This 
kind of approach is misguided as both (1) and (9) should be used in conjunction. 
From a more philosophical perspective we might say that the wavefunction VU 
belongs to the realm of quantum potentialities and is dictated by (1). On the 
other hand its absolute square which is directly measurable and seems to lay in 
a higher plane of ontological density is governed by (9). Yet, it is the wavefunc- 
tion evolution equation that guides the probability density through the velocity 
field. In the de Broglie-Bohm picture the Schrödinger equation guides the unob- 
servable quantum trajectories through the same velocity field and here we have 
moved one step further. Naturally, a question arises about the solutions that 
(9) admits. It is evident that n = |W|? satisfies the quantum Lighthill equation 
but this does not necessarily exhaust every conceivable solution, it is possible 
that there are distinct solutions of the latter such that n 4 |W|?. In that case we 
would have deviations from Born’s rule and the quantum equilibrium hypoth- 
esis would not hold. This possibility has been argued by various authors [19] 
using different kinds of argumentation from the one presented here. It should 
be remembered nevertheless that not every solution of an equation of motion, 
even an exact one, must inevitably be instantiated in nature since it could be 
unstable [20] and sensitive even to slight disturbances to remain observable. 
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